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Abstract 

A Response Function Theory and Scattering Theory applicable to the study of physical properties 
of systems driven arbitrarily away from equilibrium, specialized for dealing with ultrafast processes 
and in conditions of space resolution (including nanometric scale), are presented. The derivation 
is done in the framework of a Gibbs-style Nonequilibrium Statistical Ensemble Formalism. It 
is shown the connection of the observable properties with time and space-dependent correlation 
functions out of equilibrium. A generalized fluctuation-dissipation theorem, which relates these 
correlation functions with generalized susceptibilities is derived. It is also presented the method, 
useful for calculations, of nonequilibrium-thermodynamic Green functions. A couple of illustration 
with application of the formalism, consisting of the study of optical responses in ultrafast laser 
spectroscopy and Raman Scattering of electrons in III-N semiconductors (of "blue diodes") driven 
away from equilibrium by action of electric fields of moderate to high intensities, are described. 
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I. INTRODUCTION 



The renowned Ryogo Kubo once stated that " statistical mechanics has been considered a 
theoretical endeavor. However, statistical mechanics exists for the sake of the real world, not 
for fictions. Further progress can only be hoped by closed cooperation with experiment" [1]. 
This is nowadays particularly relevant because the notable development of all the modern 
technology, fundamental for the progress and well being of the world society, poses a great 
deal of stress in the realm of basic Physics, more precisely on Thermo-Statistics. Thus, on the 
one hand, we do face situations in electronics and optoelectronics involving physical-chemical 
systems far-removed-from equilibrium, where ultrafast (pico- and femto-second scale) and 
non-linear processes are present. Further, we need to be aware of the rapid unfolding of 
nano-technologies and use of low- dimensional systems (e.g., nanometric quantum wells and 
quantum dots in semiconductors heterostructures) [2]. All together this demands having an 
access to a statistical mechanics being efficient to deal with such requirements. On the other 
hand, one needs to face the study of soft matter and fluids with complex structures (usually 
of the average self-affine fractal-like type) [3]. This is relevant for technological improvement 
in industries like, for example, that of polymers, petroleum, cosmetics, food, electronics and 
photonics (conducting polymers and glasses), in medical engineering, etc. Moreover, in 
both type of situations above mentioned there often appear difficulties of description and 
objectivity (existence of so-called "hidden constraints"), which impair the proper application 
of the conventional ensemble approach used in the general, logically and physically sound, 
and well established Boltzmann-Gibbs statistics. A tentative to partially overcome such 
difficulties consists into resorting to non-conventional approaches [4-7]. 

Since, as noticed, a most relevant objective of any nonequilibrium statistical theory is to 
provide a comprehension of the underlying physics related to the relaxation phenomena that 
can be evidenced in experiments, it needs be coupled with a response function theory. This 
is the subject of this paper, where we specifically resort to the use of a Non-Equilibrium 
Statistical Ensemble Formalism (NESEF for short) [8-11]). 

It can be noticed that nowadays two approaches appear to be the most favorable for 
providing very satisfactory methods to deal with systems within an ample scope of nonequi- 
librium conditions. They are, on the one hand, Numerical Simulation Methods [12], or 
Computational Physics. In particular, to it belongs Non-Equilibrium Molecular-dynamics 
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NMD [13], a computational method created for modeling physical systems at the microscopic 
level, being a good technique to study the molecular behavior of several physical processes. 
On the other hand, we do have the kinetic theory based on the far-reaching generalization 
of Gibbs' ensemble formalism, the NESEF [11,14]. 

The present structure of the formalism consists in an extension and generalization of 
earlier pioneering approaches, among which we can pinpoint the works of Kirkwood [15], 
Green [16], Mori-Oppenheim-Ross [17], Mori [18] and Zwanzig [19]. NESEF has been ap- 
proached from different points of view: some are based on heuristic arguments, others on 
projection-operator techniques (the former following Kirkwood and Green and the latter 
following Zwanzig and Mori). 

The formalism has been systematized and largely improved by the Russian School of 
statistical physics, which can be considered to have been initiated by the renowned Nicolai 
Nicolaievich Bogoliubov (e.g., see ref. [20]) and we may also name Nicolai Sergeievich Krylov 
[21], and more recently mainly through the relevant contributions of Dimitrii Zubarev [8,9], 
Sergei Peletminskii [22], and others. We present in Refs. [11] a systematization, as well as 
generalizations and conceptual discussions, of the matter. 

It may be noticed that these different approaches to NESEF can be brought together 
under a unique variational principle. This has been originally done by Zubarev and Kalash- 
nikov [23] and later on reconsidered in Refs. [9,11]. It consists on the maximization, in the 
context of information Theory, of Gibbs statistical entropy (that is, the average of minus 
the logarithm of the statistical distribution function [24,25], which in Communication The- 
ory is Shannon informational entropy [26,27], subjected to certain constraints and including 
non-locality in space, retro-effects, and irreversibility on the macroscopic level. 

Concerning Response Function Theory, the usual theory to calculate linear responses 
to mechanical perturbations (e.g. [28-33]) is based on expansions in terms of correlation 
functions in equilibrium. As initial condition is taken that of equilibrium with a thermal 
reservoir, and next it is studied the evolution of the system as if it were isolated from 
all external influences except the driving field. Let us now consider the situation when a 
mechanical perturbation is applied on an already far-from-equilibrium system, in which are 
unfolding irreversible processes which are describable in terms of equations of evolution for 
a basic set of macrovariables in the non-equilibrium thermodynamic space of states. Since 
NESEF provides a seemingly powerful method to obtain a description of the macrostate 
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of such systems, it is appealing to derive a response function theory based on correlation 
functions in the unperturbed nonequilibrium state of the system. Schemes of this type have 
been proposed [31-33], and next we systematize and extend this treatment, in such a way to 
allow the treatment of experiments involving time-resolution (including the ultrafast time 
scale of pico- and femto-seconds), and space resolution (including those in the emerging 
nano-science and technology). 

It is shown the connection of the observable properties with correlation functions out of 
equilibrium; a generalized fluctuation-dissipation theorem - relating correlation functions 
and generalized susceptibilities - is derived, and the method, useful for calculations, of 
nonequilibrium-thermodynamic Green functions is presented. This is done in Sections II, 
III and IV, and in Section V we present a scattering theory, in the same conditions, namely, 
including time and space resolution, for systems far-from equilibrium. The connection of 
the scattering theory with response function theory follows from application of the nonequi- 
librium fluctuation-dissipation theorem. 

Finally, in Section VI we present a couple of illustrations showing the working of the 
theory in the study of two kind of experiments, namely optical responses in ultrafast laser 
spectroscopy of polar semiconductors, and Raman scattering of electrons in doped III-N 
semiconductors ("blue diodes") in the presence of electric fields with moderate to high 
intensities. In the latter case the nonequilibrium fluctuation-dissipation theorem allows to 
connect the Raman spectrum with nonlinear transport properties in these materials (nonlin- 
ear and time-dependent conductivity and diffusion coefficient, and a generalized - nonlinear 
and time-dependent Einstein- relation). 

II. RESPONSE FUNCTION THEORY FOR FAR-FROM-EQUILIBRIUM SYS- 
TEMS 

We consider an open many-body system out of equilibrium, which is in contact with a 
set of reservoirs and under the action of pumping sources. We are essentially presenting the 
most general experiment one can think of, namely a sample (the open system of interest 
composed of very-many degrees of freedom) subjected to given experimental conditions, as 
it is diagrammatically described in Fig. 1. 

In Fig. 1, the sample is composed of a number of subsystems, <x,-, (or better to say 
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FIG. 1: Diagrammatic description of a typical pump-probe experiment in an open dissipative 
system. 

subdegrees of freedom, for example, in solid state matter those associated to electrons, 
lattice vibrations, excitons, impurity states, collective excitations as plasmons, magnons, 
etc., hybrid excitations as polarons, polaritons, plasmaritons and so on). They interact 
among themselves via interaction potentials producing exchange at certain rates, r^, of 
energy and momentum. Pumping sources act on the different subsystems of the sample - 
via particular types of fields, electric, magnetic, electromagnetic, etc. - well characterized 
when setting up the experiment, and there follows relaxation of the energy in excess of 
equilibrium the system is receiving to the external reservoirs, Tjr. Finally, the experiment 
is performed coupling an external probing source, characterized in the figure by P(t), with 
one or more subsystems of the sample, and some kind of response, say R(t), is detected by 
a measuring apparatus (e.g. ammeter, spectrometer, etc.). 

It needs be understood that the pumping sources exert their influence on the given open 
system through the fields they generate, say, magnetic, electric, electromagnetic as produced 
for example from a laser machine, and so on, or, eventually, in scattering experiments is the 
interaction potential with the particles of an incoming beam. 
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Let the Hamiltonian be 

H(t) = H a + V(f) , (1) 

that is, the Hamiltonian H a of the system in the presence of the fields of the pumping sources 
which drive it away from equilibrium, plus the interactions with the reservoirs, and with the 
field(s), V(t), created by the external perturbing apparatus. We take for the latter the form 

V(*) = "/ d 3 r'T(r',t)A(v'), (2) 

where J 7 is the expression for the perturbing force (J 7 = —SV/SA(r'), where S is a functional 
derivative) and A(r') an observable of the system to which it is coupled. We recall that the 
system is in contact with ideal reservoirs, and the statistical operator is taken as a product 
of the one of the system g £ (t) times the stationary canonical distribution of the reservoirs qr, 
which we write TZ e (t) = g e (t) x g R [11] (see Appendix A). Moreover, H a = H + Hi + W = 
Hq + H', introducing H' — Hi + W, where Hi accounts for all the interactions in the system 
and W for the interaction with the surroundings. 
Schrodinger equation for this system, i.e. 

ih^W)) = H(t)W)) , (3) 

with the initial condition \ip{ti)) at time ti when the perturbation is switched on, has the 
formal solution 

\m) = U(t,UM(U)) , (4) 
where U is the evolution operator satisfying that 

ih^U(t,U) = H(t)U(t,ti) , (5) 

with U(ti,ti) = 1 (the unit operator), and we recall that it is a unitary operator, that is 
WU = 1. 

Let us now introduce the interaction representation, writing 

U(t,t i ) = U tr (t,t i )U'(t,t i ) , (6) 

where 

' o cwi^t-U)^} . (7) 
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which is the evolution operator in the absence of the perturbation potential V; this means 
that we are separating the internal dynamics of the system from the dynamical effects of 
the perturbation which are accounted for by U' . Introducing Eq. (6) in Eq. (5) and using 
Eq. (7), it follows that U' satisfies the equation 

ih^-U'(t,t i )=V{t)U'(t,t i ) , (8) 

with U'(ti,ti) = 1, and 

V(t) = Ul(t,t t )V(t)U a (t,t t ) . (9) 
Equation (8) can be transformed in an equivalent integral equation, namely 

- 1 t* ~ 

U'(t,U) = 1 + -J dt'V(t')U'(t\U) , (10) 



which admits the iterated solution 



U'(t,U) = 1 + 4 f dt '^( t ') + 7^ f dt ' f dt"V(t')V(t") + 
? " Ju Jti Ju 



n-1 



= 1 + £t^ dh--. dt n V{t 1 )---V{t n ). (11) 

n=1 Ju Ju 

The quantum-mechanical expected value at time t of the observable A(r), to which the 
external field is coupled, is then 

a(r,t) = m)\MrM(t)) 

= ^{u)\u'\t, tjuKt, tJAtfU.it, Wit, U) \il>(U)) 

= m^\t,t i )A{v,t)U\t^(ti)) , (12) 

where 

A(r,t) = Ui(t,t i )A(r)U a (t,U) . (13) 
According to Eq. (11) we have that 

a(r,t) = (^)|(i-4 / dt'V(t') + ---)A(r,t) x 

lh Ju 

(i + 1 f t dt'V(t') + ---)m i )) , (14) 
lh Ju 

where we have used that V is Hermitian. 



7 



Considering a weak perturbation - characterized by V(t) - imposed on the initially (at 
time U) far-from-equilibrium system, we truncate the series of terms in Eq. (14) in first 
order in V, that is, we consider from now on a linear response theory, to obtain that 

a(r,t) = a(r,U) + 1 f dt' (iP(U)\[A(r,t), V(t')MU)) , (15) 

where 

a(v,U) = m^AivMiU)) (16) 

is the expected value of the observable at time tj prior to the application of the perturbation. 
Using Eq. (2) we have that 

Aa(r,0 = -1 jfrf^y d\'mMMr,t),A(r'Xmr',t'mu)) , (17) 

where Aa(r, t) = a(r, t) — a(r, tj) is the departure of the observable A from its value at the 
initial time when under the action of the perturbing potential. 

Introducing the statistical operator for the pure (quantum mechanical) state, namely 

v(t i ) = m i ))&(t i )\, (is) 



which, we recall, is a projection operator over the vector state \ip(ti)), we can rewrite Eq. 
(17) as 

Aa(v,t) = —J* dt' J d 3 r'Tr{[A(r,t),A(r',t')]V(t t )}J(v',t') , (19) 

and where we have considered adiabatic application of the perturbation taken as U — > — oo. 

Next step is going over the macroscopic state taking the average over the nonequilibrium 
ensemble of pure states, compatible with the macroscopic conditions of preparation of the 
sample. In the usual way, if we call V n (ti) the statistical operator for the pure state in, say, 
the n-th replica, and p n the probability of such replica in the corresponding Gibbs ensemble, 
the statistical average over the ensemble of mixed states of the system, and the average over 
the states of the reservoir (system and reservoir are coupled via the interaction W) is 



A(i(r)|i> = -1 f dt' fd*r'J2PnX 

Tr{[A(r,t),A(r',t')]V n (t t ) x g R }T(r',t') 
= -j^y dt' J d 3 r'Tr{[A(r,t),A(r',t')]g £ (t t ) x g R }T(r',t') , (20) 
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with 

Qe {U) = Y.Pn'PniU) = ^2pn\MU))(MU)\ , (21) 
n n 

and we have introduced 

A(A(r)\t) =Tr{A(v)n £ (t)}-Tr{A(v)n £ (U)} , (22) 

where, we recall, lZ £ {t) = g £ {t) x g R . 
But Eq. (20) can be rewritten as 

A(i(r)|i> = -1 jf* dt' J d 3 r,Tr{[A(v),A(r',t'-t)]n £ (t)}j(v',t') , (23) 



since 



Tr{[C/t(i,^i(r)C^(t,*0,^(t , ,ti)i(r')C/ ff (t , ,* i )]7e e (* i )} = 

Tr{A(r),A(r',t'-t)K £ (t)}, (24) 

where we have used the invariance of the trace operation by cyclical permutations and the 
group properties of operators U . We stress that in nonequilibrium conditions there is no 
invariance for translation in time as a result of the time-dependence of the statistical operator 
which describes the irreversible evolution of the macroscopic state; the operator A(t' — t) in 
Eq. (24) shows time-translational invariance because such dependence on time arises out of 
the microdynamical evolution governed by Hamiltonian H a . 
Moreover, introducing the definitions 

X "(r,v';t'-t\t) = ^Tr{[A(r),A(v',t'-t)]n £ (t)} , (25) 

X (r, r'; t - t'\t) = 2i0(t - t') X "(v, r'; t' - t\t) , (26) 
with 9 being Heaviside's step function, we have that Eq. (20) becomes 

dt' d 3 r'x(r,r'-t-t'\t)T(r',t'). (27) 

-oo J 

Taking into account the integral representation of Heaviside step function, namely 

9(t) = lim % \ , 28 

we find, after some calculus, that the Fourier transform in time r = t — t' of the nonequilib- 
rium generalized susceptibility of Eq. (26) is given by 

X r, r'; co t) = / — A } \ ' , 29 

J_ QO IX to — OJ — IS 



where, we recall, s — > 40, and x"( r ) Y '\ u'\t) is the Fourier transform at frequency u of the 
r-dependent x" of Eq. (25), namely 



oo 



X "(r,r>|i) = / dr X "(r,r'; -r\t)e^ . (30) 

J — oo 

Using the fact that 

}%^h-s =pv l Tt7r6{xh (31) 

which are the so-called advanced and retarded Heisenberg delta functions and pv stands for 
principal value, Eq. (29) becomes 

X(r,r';u\t) = Re{x(r, r'; u\t)} + ilm{x(r, r'; u\t)} 

= pv + %x (r,r ;w t , (32) 

J-oo 7T U/ - W 

where Re and Jm stand for real and imaginary parts respectively, and then 

Imx(r,r';u\t) = x"(r,r' ] u\t) , (33) 

Rex{r,r;uj\t) =pv , (34) 

J-oo 71 W -u 

after recalling that x"( r ) r 'i w l^) is a re& l quantity as shown in Eq. (51). 

This Eq. (34) is one of the so-called generalized Kramer- Kronig relations, with the other 
being 

r ( , r ^ fie X (r,r>^) 

imx(r,r ;u\t) = -pv , (35) 

J-oo 7T u/ - u; 

obtained from Eq. (34) once it is used the operational relation 

/oo J // 
— (w - a/TV - a/)" 1 = -vr5(w - w') . (36) 
-oo 7T 

We recall that Kramer-Kronig relations are a consequence of the principle of causality, and 
involving the fact that x, once extended to the complex ^-plane (z = ou + iy), i.e. x( r > r 'i z \t)-> 
has poles in the lower z-plane and it is regular in the upper z-plane. Furthermore we notice 
that Rex is an even function of uj while Imx is an odd one. Thus, x( r ! r '; ^\t) nas the same 
properties as the equilibrium one presented in Eq. (52) below, as it should. 

Moreover, 7mx(r, r'; u\t) is related to the power absorption by the system. First, we no- 
tice that the external force applied on the system can be Fourier analyzed in time, obtaining 
a linear superposition of Fourier components, namely 

WO = / ^W u)(e*-* 4 e-***) , (37) 
10 



so it is a real quantity, and let us calculate the average over time t' (and then r) of the 
quantity which represents the power absorbed by the system, namely the average over a 
time interval At (typically the experimental resolution time) given by 

rt+At JeV+A 1 rt+At 



But, we do have that 

Tr{H{t')^n £ {t')) = ±Tr{[H(t>),H(t')]K E (t')} = , (39) 



dt' y ') ih 
where we have used that 1Z £ satisfies Liouville equation. Hence 



= ~ df J d\'Tr{A{r')n £ {t')} d -^^ , (40) 

with V given in Eq. (2). Taking into account that 

Tr{A(r')n £ (t)} = Tr{A{r')U{t\U)n £ {U)U\t\U)} 
= Tr{U ] (t', ti)A(r')U(t', t^K^U)} 

= Tr{U'\t',U)A(r',t')U'(t',U)n £ (U)} , (41) 

where we used Eqs. (6) and (13), and next, resorting to Eq. (11) in first order (linear 
response) and to Eq. (20), we find that 

rt+At 



W(t) = ^ jf dt' J d 3 r'j dt" J d 3 r"Tr{[A(r',t'),A(v",t")} 



x 



n £ {u)}H^"^ d:F{j ' )tr) 



dt 

"t+At 



= ~hl dt 'S rfV ^J)^r- ' (42) 

with A A of Eq. (20). Because of Eq. (27) we can write Eq. (42) in the alternative form 

r-t+At 



W(t) = ~J dt' J d 3 r' y dh" y dr X (r>, r"; r\t') 



x 



Wr+0 «?m, (43) 

n 



with t = t" — t 1 ', and using Eq. (37) and the Fourier transform of x, it follows 



1 rt+At 



{x-*(r',r"; W |(')(e 2 '"'' - 1) - X '(r', r"; (e" 2 '"'' - 1)} . (44) 

If we consider a particular situation in which lZ E (t') varies weakly in time in the interval 
At (implying in that ultrafast relaxation processes are not present), we can approximate it by 
TZ £ (t), and in conditions such that At involves several periods 27r/w so that the exponentials 
cancel on average, using Eq. (23) it follows that 

W(t) = X - J d 3 r' J d 3 r" J duT(r', u)T(r", u)a(r', r"; u\t) , (45) 

where 

a(r',r";u\t)=u X "(r',r";u\t) (46) 

can be considered a kind of absorption coefficient at each frequency uj associated to the per- 
turbing force, and at the macroscopic (nonequilibrium thermodynamic) state of the system 
at time t. 

Evidently, and this is a fundamental point to be stressed, the generalized susceptibility 
depends on g £ , through TZ £ , and since the latter is a functional of the time- (and eventually of 
the space-) dependent variables that characterized the nonequilibrium thermodynamic state 
of the system, then the calculus of responses needs be coupled with the one of the equations 
of evolution for the basic variables, as given by the nonlinear quantum kinetic theory that 
the formalism provides [8-11,14]. 

Closing this section let us add some additional considerations. 

First, if the system has translational invariance (or near translational invariance as in 
the case of regular crystalline matter), the dependence of x" an d X on r an d r' is through 
the difference r — r'. We can then introduce the Fourier transform in space, namely 

X"(k,u,|t) = J d 3 b X "(bMt)e- ikh , (47) 

where b = r — r'. Moreover, 

[x"(r,r>;-T\t)r = -x"(r,r';-T\t), (48) 

as a result of x" involving a commutator of Hermitian operators, and this x" is then a purely 
imaginary quantity. Similarly, it follows that 

X "(ry;-r\t) = -x"(r,r';r\t), (49) 
12 



what implies that 

X "(Kco\t) = -x"(-k,-co\t) . (50) 
On the other hand, on the basis of Eqs. (32) and (50) we have that 

Im X (k,u\t) = l( X (k,w|t)-x*(k,w|t)) 

i f°° r , r ,,, / e -i(k b- w 'r) e »(k.b-«'T) 



2i ./_oo vr 



rf 3 6 y drx"(b,-r|*) 
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= x"(Koo\t), (51) 

which is then a real quantity, and also is x"( r > r '; M*)- 

Moreover, if the initial condition of preparation of the system is the one of equilibrium 
with the reservoirs, characterized by distribution g cq (which commutes with H), we recover 
the usual expression [28] 

X (r, r'; t - t'U = 6(t - t') l -Tr{[A(v), A(r', t' - t)] Qeq } . (52) 

where, as usual, because the equilibrium has been established, the interaction of system and 
reservoir is neglected. 

Secondly, as pointed out in the Introduction the notable developments in instrumentation 
that are at present being accumulated, and which are necessary for the study of systems 
working in far-from-nonequilibrium conditions in the sought-after miniaturized devices with 
ultrafast responses of nowadays advanced technology, require the mechanical-statistical anal- 
ysis in short time intervals (as described above) and in nanometric spatial regions. This is 
also contained in the theoretical treatment already described. In fact, if the property of 
the systems is measured in a small region around position r, and also evolving in time, the 
expected value is given in Eq. (27), which we can alternatively write as 

A(A(r)|i> = -J Jt'Tr{[A(v,t),V(t',t' -t)]n E (t)}} , 

with [cf. Eq. (2)] 

V(t', t'-t) = - J d 3 r'F(r', t')A(r', t' - t) , 

Finally this NESEF-based approach to Response Function Theory involves, as seen above, 
the calculation of averages in terms of the nonequilibrium statistical operator, a quite difficult 
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task. First we recall that, we have taken the system in contact with external reservoirs, 
the latter very much larger than the system, which, for all practical purposes, remain in a 
stationary state of equilibrium all along the realization of the experiments, and that lZ e (t) = 
6e{t) x Qri where qr is the stationary equilibrium statistical operator of the reservoir(s), and 
g £ {t) the nonequilibrium statistical operator of the system. Moreover, we can write [8-11] 



that is, the sum of the "instantaneously-frozen" auxiliary statistical operator g and the con- 
tribution g' £ which accounts for the relaxation processes developing in the media. Therefore, 
we can write Eq. (25) as 



where in x" the averaging is over the auxiliary ensemble, characterized by g, and x" is the 
averaging in terms of the contribution g' £ . Using the expression for g' £ , given in terms of 
g [11], it can be shown that it is expressed in a Born-perturbation-like series in powers of 
H', the internal interactions in the system. Therefore, whereas the weak-coupling limit can 
be used, we can retain only x to a good degree of approximation. In the NESEF-based 
kinetic equations the use of such limit renders the equation Markovian in character [14]. 
We stress again that the Response Function Theory for systems away from equilibrium is 
always coupled with the kinetic equations that describe the evolution of the nonequilibrium 
thermodynamic state of the system. 

Let us next see another important property of the nonequilibrium generalized suscepti- 
bility, namely a fluctuation- dissipation theorem in far- from- equilibrium conditions. 

III. FLUCTUATION-DISSIPATION THEOREM IN FAR-FROM-EQUILIBRIUM 
CONDITIONS 

The fluctuation-dissipation theorem (fdt) - originally a relation between the equilibrium 
fluctuations in a system and the dissipative response induced by external forces - provided 
a major impetus for the development of discussions of irreversible processes. A classical 
particular form seems to have been proved by Nyquist [34] for the relationship between the 
thermal noise and the impedance of a resistor. Derivation from phenomenological points 



g £ (t) = g(t) + g' e (t) , 



(53) 



X "(r, r'; t' -t\t) = x"(r, r'; t' -t\t) + ^'(r, r'; t' -t\t). 



(54) 
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of view followed, together with stochastic approaches, and finally entered into the domain 
of statistical mechanics [35,36]. We here extend Kubo's approach in order to encompass 
arbitrary nonequilibrium conditions. 

As noticed those results involve the immediate neighborhood of the equilibrium, and 
were very well established. For systems out of equilibrium (particularly those far from 
equilibrium) the situation is not clearly delineated, and some approaches are available for 
steady-state conditions in a stochastic approach [37] and in transient regimes for particular 
ensembles [38] . We address here the derivation of a fdt for far-from-equilibrium systems, in 
the framework of the nonequilibrium ensemble formalism NESEF, which is a generalization 
to arbitrary nonequilibrium conditions of the formalism developed by Kubo [36] in the case 
of systems in equilibrium. 

A fluctuation-dissipation theorem for systems arbitrarily away from equilibrium, and in 
the formalism of NESEF here presented, follows from the comparison of two expressions: one 
is a correlation function of two quantities and the other a dynamic response of the system to 
an external deterministic perturbation, that is, the generalized susceptibility of the previous 
section. 

Let us first recall the case of equilibrium [36]. Consider the quantities A and B: their 
correlation function over the canonical ensemble in equilibrium is given by 

S AB (r, r', t-t')= Tr{AA(r, t)AB(r', t')g c } = Tr{AA(r, t - t')AB(r')g c } , (55) 

AA = A — Tr{Ag c }, etc., and the generalized susceptibility is 

X 'Ur,r>;t-t>) = ±_Tr{[A(r,t), B(r',t')] Qc } = ±_Tr{[A(r,t - t f ), B(r')] Qc } , (56) 

which is a generalization of the one of Eq. (25), with 

A(T,t) = e-fc 6 Ae* t6 , (57) 

where H is the Hamiltonian of the system (in the absence of any external perturbation), 
that is, the operators are given in Heisenberg representation, and 

H 

e k s T 

Qc = Z(T,N,V) (58) 
is the canonical distribution in equilibrium. Using the operational relationship 

e-^ H e-rn^A(r)er^ = e ~^( W A{r)e^ t+i ^ A = A(r, t + i(5h)e~^ , (59) 

15 



we do have that 

Tt{q c A(t, t)B(v\ t')} = Tr{g c B(r', t')A(r, t + ih/3)} . (60) 

Furthermore, let us introduce the Fourier transform in time of the correlation function, 
namely 

S AB (r, r'; t - f) = f°° pS AB (v, r'; co)e^^ , (61) 

and consider 



/oo 
rfrTr{A J B(r / )Ai(r,r) f?c } ( 
-oo 

/oo 
dTTr{£ c A£(r')Ai(r,T + i/m)} , (62) 
-oo 

where we have used the relationship of Eq. (59), and r = t — t'. Introducing r = r' + i/3h it 
results that 

S AB (r', r, -u) = S BA (r, r'; w ) e -^ . (63) 

On the other hand 

r°° i 

/•°° 1 
= y rfre^-Tr{[AA(r,r),A J B(r')]&} 

= ^{^(r, r'; u) - S BA (r', r; -a;)} , (64) 

where we can introduce AA and AB because the extra terms cancel in the commutation, 
and using Eq. (63) it follows that 

S AB (r, r'; u ) = 2h[l - e'^x'Ur, r'; u) , (65) 

which is the traditional form of the fluctuation-dissipation theorem. In this condition in the 
linear regime around equilibrium, the dependence on the space coordinates is usually of the 
form r — r', and then making the Fourier transform in the space variable we have that 

S AB (k,u) = 2h[l - e-^rV; B (k,u;) , (66) 

where x" is the imaginary part of the generalized susceptibility of the previous Section [cf. 
Eq. (33)]. 
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Let us now go over the case of a system away from equilibrium, defining 

S AB (r,t;r',t'\t l )=Tr{AA(r,t)AB(r',t')g £ (U)} , (67) 

X %(v,t;r',t%) = ^Tr{[A(r,t),B(v',t')]g £ (U)} , (68) 

where g £ (ti) is the distribution characterizing the preparation, in nonequilibrium condi- 
tions, of the system at time U when the experiment is initiated, and AA(r, t) = A(r, t) — 
Tr{A(r,t)g £ (ti)} = A(r,t) - Tr{A(r)g £ (t)}, after using that 

A(r,t) = rf(t)A(r)U(t), (69) 

Q e (t) = U(t)Q e (ti)tf(t), (70) 
U(t) = e-fo*-^ 6 . (71) 

Further, we write for g £ 

g £ (t) = exp{-Se(t)} , (72) 

where 

S e (t) = S(t,0) + Ut), (73) 

with 

Ut) = - f dt> ±S(lf, t'-t). (74) 

Proceeding along a similar way as done in the case of equilibrium, we first take into 
account that we can write 

Tr{AS(r,t)AA(r', *')&(*<)} = Tr{U\t)AB(r' ,t)U(t)U\t)AA( r )U(t)g £ (t t )} 

= Tr{AB(r',t)AA(r)g £ (t)} . (75) 
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Moreover 

Tr{A(r,t)B(r',t > )g £ (t t )} = Tr{B(r' ,t')g £ (ti)A(r,t)} 

= Tr{e-^(*^e^ ( * !) f/ t (t , )5(r')f/(t')e'^ ( ' l) f/ t (t)i(r)f/(t)} 
= Tr{e~^ u) U\t)U(t)e^ {u) U\t)U{t)U\t') x 
B(v')U\t')U\t)U{t)e-^ {u) U\t)A{v)U{t)} 
= Tr{(U(t)e- S <Wu\t))(U(t)e S <teW(t)) x 



{U{t)U\t')B{v')U{t')U^){U{t)e- s ^U\t))A{r) 

= Tr{e- § ^e § ^B(r',t' -t)e- § ^A(r)} 
= Tr{B £ (r',t'-t\t)A(r)g £ (t)}, (76) 



where 



4(r', -r\t) = e^ (t) U(T)B(r')U\r)e-^ (t) , (77) 

with t — t — t\ and we recall that U(—t) = U^(t). 
Consider now the susceptibility, when we have that 

X^M;r',^) = ^Tr{[A(r,t),B(r / ,t')]g £ (t t )} 

= ^Tr{[Ai(r,t),AB(^,0]&(*i)} 

= ^Tr{(Ai(r,t)Afi(r',f) - Afi(r',f)Ai(r,t)) & (t i )} 

= ^Tr{ (A4(r', -r\t)AA(r) - AB(v' , -r)AA(v))g £ (t)} , (78) 

where we have used Eq. (76), and it can be noticed that the distribution g £ in the last 
expression can be given as given at time t when a measurement is performed. Equation (78) 
can now be written in the form 

X %(r, r'; t' - t\t) = ^ [S BA (v', r; t' - t\t) - S BA (r', r; t' - t\t)\ , (79) 

after defining 

S BA (r', r; t> -t\t)= Tr{AB £ (v f , -r) AA(v) g £ (t)} , (80) 

S BA (r', r; t' - t|t) = Tr{A£?(r', -r)AA(r)g £ (t)} . (81) 

Taking into account that the nonequilibrium distribution admits a separation consisting 
of the addition of two parts, one is the so-called relevant part, g, plus a contribution, g' £ , 
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which accounts for relaxation processes which are governed by H' in Eq. (1), we separate 
the fluctuation-dissipation relation of Eq. (64) in a "relevant" part (the one depending on 
q and H alone) and the rest. That "relevant" part is then 

xUr,r\-r\t) = ±Tr{ (e^ AS(r', -r)„ e -^Ai(r) - 

AB(r',r) o AA(r))g(t,0)}, (82) 

where 

AB(r, -r) = e^ T *°B(r)e-^ TA ° . (83) 

Transforming Fourier in r and also in the space coordinates after assuming dependence 
on r — r', that is, in the cases when the system displays translational invariance, we have 
that Eq. (79) becomes 

X 'UKoo\t) = ±[S% A 0t,u>\t) -S BA (k,co\t)} , (84) 

constituting a fluctuation-dissipation relation for systems arbitrarily deviated from equilib- 
rium. In particular, it goes over to the one in equilibrium [cf. Eq. (66)] when g £ (t) is 
substituted by the canonical distribution g c , as it should. 

On the other hand, in the case of space- and time resolved experiments in systems with no 
translational invariance then for the response in region Ar (the experimental space resolution 
with, say, micrometer or nanometric scale) around position r, we do have 

XAm (r;t'-t\t) = Jd 3 r'x' AB (r,r';t'-t\t) = 

= ^\S e BM {v;t' -t\t) - S BM {v;t' -t\t)} , (85) 

(with I for local), with the last two correlations over the nonequilibrium ensemble corre- 
sponding to the integration over r' of those on the right of Eq. (79). 

We illustrate these results for a specific case, say, a gas of free fermions for which 

H = ^2 e k4 C k » (86) 

k 

where c(c^) are the annihilation (creation) operators in state k (we omit the spin index), and 
ek is the energy- dispersion relation. We consider the case of the nonequilibrium generalized 
grand-canonical ensemble [11], when the set of basic variables present in the nonequilibrium 
statistical operator are the densities of energy and of particles, the flux of both (currents), 
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and all the other higher-order tensorial fluxes of order r = 2,3,..., with r being also the 
tensor rank. Further, we restrict the analysis to the homogeneous condition, i.e. densities 
and fluxes do not depend on the space coordinate, and then, in reciprocal space, 

S(t, 0) = W F h (t)e k + F n (t) + ^T[if ] (t) ® u M (k)e k + F^{t) ® «W( k )] b Ck ( 87 ) 

k I r>l J 

is the informational entropy operator in this case. 

In this Eq. (87) the F's are the intensive nonequilibrium thermodynamic variables, 
associated to the energy, the number of particles, and the fluxes of these two of order r 
(=1,2,... ); -uM(k) is the tensor consisting of the tensorial product of r times the generating 
velocity u(k) = ft^Vkek, namely the group velocity of the fermion in state k; sign <g> stands 
for fully contracted product of tensors. 

Equation (87) can be written in the compact form 

S(t,0) = ^y9k(t)c k c k , (88) 

k 

where, evidently, ip is the quantity between curly brackets in Eq. (87). Hence, in these 
conditions we have that the "relevant" part of the correlation of Eq. (80) is 

= Tr|e^ EkTkekC k Ck 5(rOe-^^ rk£kC k Ck i(r)^(t,0)| , 



(89) 



where we have defined the quantity 

Tk = T-^(y9k(t)/e k ) , (90) 

which has dimension of time. 

Let us further simplify matters and use a truncation in the basic set of macrovariables, 
retaining only the energy Hq, the number of particles N, and the flux of matter I n , which 
multiplied by the mass of the fermions becomes the linear momentum, P. We call the 
associated nonequilibrium thermodynamic variables F h = j3*(t), F n = — (3* (t) /i* (t) , and 
F„(t) = — f3*(t)m*v(t)i introducing /3* _1 (t) = k B T*(t), a reciprocal of a quasitemperature, 
yu*(t) a quasi-chemical potential, and a drift velocity v(i) [11,39-42]. This implies in that 
we are using a kind of nonequilibrium canonical distribution with an additional term arising 
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out of the presence of the current with drift velocity v(i). Moreover, let us choose for A and 
B the nondiagonal elements of Dirac-Landau-Wigner single-particle density matrix written 
in second quantization, namely 

A = 4 Q = C LiQ C k+|Q > ( 91 ) 

B = h kQ = 4 + i Q c k _i Q , (92) 

which will appear in the calculation of inelastic scattering cross sections later on. Hence the 
relevant part of the corresponding correlation is 

^ nrit (k,Q,-r|t) rcl =Tr{U(T\t)n kQ U\T\t)ht Q g(t,0)} , (93) 

where 

U(r\t) = expj^ [ Ck ( r - - ihP*(t)n*(t) + ihp*(t)v(t) ■ Hk] c[c k } , (94) 



and we have that 

U(r\t)n kQ W(r\t) = expj i(-r + i^*(*))(e k+ i q - e k _i Q ) + ih(3*(t)v(t) • Q Jn kQ 

= e /3 *W v W' Q n kQ (-r + ^/3*(t)) (95) 

Hence 

5 nnt (k, Q; -r|t) rcl = e h ^^Tr{n^{-r + thf3*(t))nl Q g(t, 0)} , (96) 
and using Eq. (82), after some calculation we find that 

5 nnt (k, Q; -r\ty d = 2h[l - e-«^- v W-Q)]- 1 Xn „ t (k, Q; co\ty cl (97) 
For v = and in the case of equilibrium we recover the well known result of Eq. (66). 



IV. NONEQUILIBRIUM-THERMODYNAMIC GREEN FUNCTIONS 

According to the previous sections, to obtain response functions requires the calcula- 
tion of nonequilibrium correlation functions. This is a difficult mathematical task which 
can be facilitated by the introduction of appropriate nonequilibrium thermodynamic Green 
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functions [43-45]. The approach is an extension of the equilibrium thermodynamic Green 
function formalism of Tyablikov and Bogoliubov [46] . 

We define the retarded and advanced nonequilibrium-thermodynamic Green functions of 
two operators A and B given in Heisenberg representation, by the expressions 

((A(v,r);B(v')\t))^ = ^e(r)Tr{[A(r,r), B(r%TZ E (t)} , (98a) 

({A(v,r);B(v')\t))^ = -Le(-r)Tr{[A(r, r), B(r%K £ (t)} , (98b) 

where t — t' — t and rj = + or rj = — stands for anticommutator or commutator of operators 
A and B. These Green functions satisfy the equations of motion 

ih^((A(r,r)-B(r')\t))^ 

(99) 

In Eq. (99), and in what follows, [A(r), B(r')] without subscript is the commutator of 
quantities A and B. Introducing the Fourier transform 

/oo 
rfrexp[^r]((^(r,r); J B(r')|t)) r ,, (100) 
-oo 

equation (99) becomes 

nuj((A(r);B(r')\uj;t)) v = ^Tr{[A(v), B(r%K £ (t)} + (([A(r), H]; B(r')\u;t)) v (101) 

A. Green Functions and the Fluctuation-Dissipation Theorem 

Next, we establish the connection of these Green functions with correlation functions. 
Consider the nonequilibrium correlation functions 

F AB {v,v'-T-t) =Tr{A{v,r)B{v')TZ £ {t)} , (102a) 

F BA (r,r';r;t)=Tr{B(r')A(r,r)n £ (t)} , (102b) 

and let \n) and E n be the eigenstates and eigenvalues of the Hamiltonian H. Defining the 
nonequilibrium spectral density functions 

J AB (r,r';uj\t) = 2tt ^(n|i(r') |m> (m\B(r') \l) (l\K e (t)\n)6(hu> - E m + E n ) , (103a) 

Iran 

K BA (r,r>|t) = 27r^(n|fi(rO|m)(m|i(r / )|/)(^ £ W|n)5(^-^ + E m ) , (103b) 

Imn 
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we obtain the relations 

/°° dui 
—J AB (r,r';uj\t)exp[-iujT], (104a) 
■oo 271 

/oo 
— K BA (r, r'; w|t) expf-iwT] , (104b) 
-oo 27r 



and 



((i(O;S(O|a;±w;t>> + + ((i(r0;S(r')|a;±< S ;t>>_ = i f^ Vr.^) ? (1Q5a) 

ft J-oo ^ W — W ± ?S 

(<i(0;S(i/)|a;±< S ;t>>_-((i(r');S(r')|a;±w;t>> + = i f , (i 05 b) 

ft > /_ 00 7T — ± ?S 

with s — >■ 40, plus sign is for the retarded and minus sign for the advanced Green functions, 
and we made use of the relation 



j — c 



dT9(±r)exp\i(uj -cu')t] = ± — — . (106) 

0J — UJ ± IS 



Equations (105) may be considered particular generalizations of the fluctuation- 
dissipation theorem for systems arbitrarily away from equilibrium. Near equilibrium, re- 
placing g £ by the canonical Gibbs distribution we recover the well known result 

{(A{v>)-B(v')\u + is))^ - ((A{v')-B{v')\u-is))^- = 1(1 - ijexpHM) , 

(107) 

where j3 = (ksT) -1 . 

We recall that the nonequilibrium-thermodynamic Green functions of Eqs. (98) depend 
on the macroscopic state of the system, and therefore their equations of motion, Eqs. (99) or 
(101), must be solved coupled to the generalized nonlinear transport equations for the basic 
set of nonequilibrium thermodynamic variables [14]. Finally, if we write for the interaction 
energy V = A exp[iut]B(r') where A is a coupling strength constant, it follows that 

(i(r)|*)-(i(r)|t)° = -1 T dreM-^r]Tr{[A{v),B{v\r)]n £ {t)} + C .c. 

in 

A f°° 

= ~ih dT9(T)[F AB (r,r';r;t) - F BA (r,r';T;t)]exp[iuT]+c.c. 
A f 00 duj'K A§ (r,r';u'\t)-J M (v,r';u'\t) | ^ 



h J _ OQ 2n ui — to' 4 is 

2\Re{{B(r'); A(r)\uj + is;t))_ , (108) 
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where Re stands for real part, and we have used the definition of the advanced Green 
function of Eq. (98). Hence, the linear response function to an external harmonic pertur- 
bation is given by an advanced nonequilibrium-thermodynamic Green function dependent 
on the macroscopic state of the system characterized by the nonequilibrium thermodynamic 
macrovariables Fj(t) [or equivalently Qj(t)], as described in Section II. 

Closing this section we note that since the nonequilibrium thermodynamic Green func- 
tions of Eqs. (98) are defined as nonequilibrium averages of dynamical quantities, recalling 
the separation of g £ in a secular and non-secular (dissipative) parts, we can write 

«i(r'); B(r>; t» = «i(r'); B(r')\u; t» scc - + ((i(r'); B(r>) \u; t))' , (109) 

where 

1 r°° . 
((A(r');B(r')\u;t)r c - = ±- / dr exp[tur]e(±r)Tr{[A(r' , r), B(r%n(t, 0)} , (110a) 

^ J —oo 

((A(r');B(r>)\u;t)y = ±\ f°° dr exp[tur]6 (±r) {[A(r' , r), B(r%; K' e (t)\t} . (110b) 

J —oo 

We notice that, in general, the last term in Eq. (101) couples the equation for the Green 
function with higher order Green functions, whose equations must be written and thus one 
obtains a hierarchy of coupled equations. Usually one solves this hierarchy of equations in- 
troducing a truncation procedure, like some kind of random phase approximation. For these 
nonequilibrium-thermodynamic Green functions a second type of expansion and truncation 
is also present, which is that associated with the irreversible processes encompassed in the 
contribution TZ' £ (t) to the statistical operator present in Eq. (110b). Care should be taken 
to perform consistently both types of truncation procedures, i.e. to maintain terms of the 
same order in the interaction strengths. We recall that the Markovian approximation in the 
NESEF-based kinetic theory [14,47] requires to keep terms containing the operator for the 
interaction energies up to second order only. The formalism of this section was applied to 
the study of time-resolved Raman spectroscopy as described in Refs. [48,49]. 

V. THEORY OF SCATTERING FOR FAR-FROM-EQUILIBRIUM SYSTEMS 

In a scattering experiment a beam of particles (e.g. photons, ions, electrons, neutrons, 
etc.) with, say, energy e and momentum Kk , incide on a sample, where they interact with 
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one or more subsystems of it (say, atoms, molecules, electrons, phonons, etc.). The particles 
are scattered, as a result of that interaction, involving a transference of energy AE, and 
of momentum hq, consequence of an excitation being created or annihilated in the system; 
Fig. 2 shows a scheme of the experiment. 




FIG. 2: Scheme of an experiment of scattering. 

Calling ei and hki the energy and momentum of the scattered particle, conservation of 
energy and momentum require that 

AE = e -e 1 , (111) 
hq = hk -hk 1 , (112) 

or 

q 2 = kl + k\ - 2k h cos 9 (113) 

after scalar product of Eq. (112) with itself (9 is the so-called scattering angle, see Fig. 2). 

Let us go over the general theory: The scattering can be characterized by the quantity 
differential scattering cross section, d 2 a(AE,q). It is defined as the ratio between the 
number of scattered particles that are collected by a detector within an element of solid 
angle dQ(9,Lp) in direction (9, </?) per unit time (to be designated SN), and the flux of 
incident particles, namely the number of particles that enter the sample per unit of time 
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and unit of area (to be designated <3> ). The latter is given by 

%=nv , (114) 

where n is the density of incident particles and vq their mean velocity (e.g. velocity of light 
in the case of photons, the thermal velocity in the case of thermalized neutrons). On the 
other hand we have that 

SN = nV w P0 ^ p[ (AE), (115) 

where w is the probability per unit time that an incident particle with momentum p makes 
a transition to a state of momentum p' l5 the latter in directions contained in the solid angle 
d£l{6, ip) whose axis we indicate by pi (or ftki); AE is the transfer of energy in the scattering 
event, and nV is the number of particles and V is the active volume of the sample, i.e. the 
region involved in the process, for example the region of focalization of the laser beam in 
the scattering of photons. 

However, since dVL is small (fixed by the size of the detector window) we can take all the 
contributions in the sum over pi as the same, i.e. to a good degree of approximation equal 
to w Po ^. Pl (AE). Hence 

d 2 a(AE,q) = ^ = -w po ^ pi (AE)b( Pl ,dQ) , (116) 

where 

b(p,dQ) = 1 - T^ p " dpidQ ( 117 ) 

is the number of states of the particles in the scattered beam, which are entering the detector, 
and we recall that p = ftk and the sum is over the plane-wave state of wavevector k. 

Therefore, for calculating the differential cross section we need to evaluate the transition 
probability per unit time, w. For that purpose let us consider a system with Hamiltonian 
H a , let us call Hp the Hamiltonian of the particles used in the experiment, and V the 
interaction potential between the system and the particles, that is 

H = H a + Hp + V . (118) 

Let us introduce the notation and |p) for the eigenf unctions of the system and the 
particles in the probe, i.e. 

H v \fi) (119) 
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H P \p) = tiu; p \p) , (120) 

where states |p) are plane waves for the free particle with momentum p in the incident and 
scattered beams. Moreover, let 

^(fO) = |*(*i))|po(*i)> (121) 

be the wavefunction at the initial time U, that is, the initial preparation of the system in 
the experiment. 

As we have seen in previous sections, it is convenient to work in the interaction represen- 
tation, and then the wavefunction of the system and probe at time t is 

|^(*)> =Uo(t,t i )U'(t,t i )\1>(t i )) (122) 

[cf. Eqs. (5) to (8)], with 

ih^U'(t,t i )=V(t)U'(t,t i ), (123) 

where V(t) is the potential in the interaction representation, i.e. evolving with H = H a +H pi 
[cf. Eqs. (9) and (10)], U is given by 

U (t,U) = UvfatjUpfaU) 

and we recall that the iterated solution of Eq. (10) is given in Eq. (11). 
Moreover, if we define the function 

m)) = ui(t,u)mt)) = u'^umu)) , (125) 

we can easily verify that it satisfies the equation 

i*^W))=V(t)W)), (126) 
with \ip{ti)) = \ip(ti)), which can be rewritten as 

and then has the iterated solution 

"t rtn-i 



\m) 



%Wr I dtl "'l dt ^(ti)---v(tr, 
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and for later convenience in the calculations we notice that we can write that 



V(t)\m) = EWMU)) , (129) 
with S, called the scattering operator, satisfying the integral equation 



S(t) = V(t) 



1 + 



(130) 



It can be noticed that the right hand side of Eq. (128) provides the effect of the pertur- 
bation V in all orders over the initial nonperturbed wavefunction, equivalent to the effect in 
first order over the interaction-representation function \ip(t)). 

Let us now fix the scattering channel, i.e. we consider, as required by Eq. (115), the 
scattering event with probe particles making transitions between states of momentum |p ) 
and |pi). According to the general theory of Quantum Mechanics, the probability for this 
event at time t is given by 

p P0 ^ P1 (*) =^i<pi, /#(*)> r, (i3i) 

where the summation over all states of the system will be a posteriori restricted by the 
selection rule involving conservation of energy and momentum in the scattering events. 
We can rewrite Eq. (131) as 

Jpo-piW = ^|(pi,/i|^o(M^'M)IV^)>l 2 

= E|( Pl ,/^ (M*M0>l 2 , (132) 

where we have used Eq. (125). On the other hand 

(p,li\U (t,ti) = (p,/i|exp|-^(t-t i )( J B M + ^; p )| , (133) 
then the exponential gives a modulus 1 in Eq. (132), and we have that 
Ppo- P (*) = £|<Pi, /#(*)> | 2 

= 5^i< Pl ,>a|[i + ^ f dt mum 2 



EKPi'/^ f dt'S(t')mu),p )\\ (134) 



i r* 

^ Kpi»a*I 
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where we have used Eqs. (127) and (129), together with Eq. (121), and the fact that the 
plane wave states |po) and |pi) are orthonormal states. 



Introducing 



S q (f) = (pi|S(f)|Po> = J d 3 r"-^S(t)—j=- , (135) 

with, we recall, hq = p — Pi, and using that the squared modulus can be written as the 
product of the complex number times its complex conjugate, we can write 

fpo-PiW = f dt"Sl(t")\n)(n\ f dt' S q (f)|<&(*i)> 

J -CO J —CO 

= W f dt ' f ^^{^1(0^(0^0^)} , (136) 

J — OO «/ — OO 

with Vo(ti) being the projection operator (statistical operator for the pure state \Q(ti)) of 
Eq. (121), and we have considered adiabatic application of the perturbation in tj — > — oo), 
implying in that the initial ultrafast transient is ignored. 

So far we have a purely quantum-mechanical calculation, and we have an expression 
depending on the initial preparation of the system as characterized by the statistical operator 
for the pure state given above. We need next to introduce the statistical average over the 
mixed state, what is done averaging over the corresponding Gibbs ensemble of all possible 
initial pure states compatible with the thermodynamic condition of preparation of the system 
at time ti, that is, we do have have that 

(Ppo-p^t)) = i /* dt' f dt'TriSKO-^Q^U) x qr} (137) 

where g £ (ti) x g R = 1Z £ {tj) is the corresponding statistical operator, involving g £ of the 
system in interaction with the thermal bath and qr of the thermal bath, which has been 
assumed to constantly remain in equilibrium at temperature T . 
The rate of transition probability w in Eq. (116), is then 

w P( ^ Pl (AE\t) = j t (P P0 ^ Pl (t)) 



= ^ J dt'Tri&ME^QM x g R } + 

1 J dt'TriS^E^QM x QR } . (138) 

Using that 

^{^(^^^[/(^^^o^'^O^o^,^)^^)} = Tr{S^(lf-t)Q e (t) x g R } , (139) 
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and that 



S^t'-t) = ( P i|4(f -t)Ul(t' -t)E^U a (lf -t)U P (t' -t)\p ) 

= e-inif-m^Po-^s^t' -t) a , (140) 

where 5q(t' — t) a = U*(t' — t, ti)EqU a (t' — t, U), we can write 

= 1 J dt' eW-QTri&E^t - t) a g £ (t) x g R } + c.c. , (141) 

where oj = ou po — uj pi , and it can be noticed that the statistical operator is given at time t, 
when a measurement is performed. 

If we consider the case of equilibrium, i.e. we take the canonical distribution g c instead 
of g e (t), and taking into account that g c and H a commute, we find that 

1 f°° 

w(q, u) cq = - I dr' e-^'TrisKr')^} , (142) 



— oo 



which is the known temperature-dependent rate of transition probability (see for example 
[50]). 

In conclusion, the differential cross section is then given by 



d 2 a(q,u\t) V 2 g(u) 1 



dudtt (2tt^) 3 v K 2 

f dt' e-W-^Tri&te -t)^g e (t) x g R } 

J -co 



f dt' e^-^TriS^H - t) a n £ (t)} + 

J — oo 

(143) 



each term within the square bracket is the complex conjugate of the other and then the 
quantity is real, as it should, and we have defined 

p 2 dp = g(u) du , (144) 

introducing the density of states g(oo) which follows in each case once it is given the dispersion 
relation ou p . 

Moreover it is stressed the fact that, differently to the case when the system is in equilib- 
rium, in the nonequilibrium initial preparation of the sample the scattering cross section is 
not closed in itself, but it needs to be coupled to the set of kinetic equations that describe the 
evolution of the out of equilibrium system, i.e., those that determine the statistical operator 
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TZ e {t). It can be noticed that this is a question also present in the response function theory 
of the previous sections. 

As it has been noticed in previous Section, we can write lZ e (t) = g £ {t) x g R and make 
use of the separation g £ (t) = g(t) + g' £ (t) to obtain that 

d 2 a(q, u\t) = d 2 a(c[, w\t) + d 2 a' £ (q,u\t) , (145) 

that is, the contribution d 2 a where the trace is taken with g and d 2 a' £ with the trace taken 
with g' £ . 

Let us consider the case of time- and space- resolved scattering, that is the detector in 
Fig. 2 collects the scattered particles arriving from an element of volume AV(r) around 
position r in the sample. For simplicity we take the first-order scattering consisting that in 
Eq. (142) we take of the scattering operator S of Eq. (130) the first contribution V, and for 
the latter we write 

N' N 

v = £ X>(r* - , ( 146 ) 

li=l j=l 

where Yj is the position of j-th particle in the system and R M the position of the /z-th particle 
in the incident beam. 
Therefore, we have that 

N 

V q = (Po|V|pi) = n b v(<l) > ( 147 ) 

3=1 

where we have introduced the Fourier amplitude 

v(q) = [ d 3 bv(b)e^ h , (148) 



with b = r — Tj, and we recall that ftq = po — pi, and rib is the density of particles in the 
beam. 

Retaining in Eq. (143) the contribution in first order in V only, it follows that 

d 2 a(q,u\t) V 2 g(u) 2 2 

dwdTl = (2^7^ Kq)l S ^Mt), (149) 

where 

S nn (q,u\t) = J2[ df e^-QTrie-^V-Q-^Qeit) x g R } + c.c. , (150) 
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We introduce now the density operator 

h(r,T) = Y,S(r-r s (t)) (151) 

j 

and then we can rewrite the correlation function of Eq. (150) as 

S nn {q,u\t) = J d 3 r J d 3 r' J dt' e^' e'^-^Trin^r, t' - t)n{r')g £ {t) x g R } + c.c. , 

(152) 

where the integrations in space run over the active volume of the sample (region of concen- 
tration of the particle beam), or in the case of a space-resolved experiment over AV^(r) and 
then we do have the time- and space- resolved spectrum 

d2 ^Zdn^ = Ay(r) / dV / ^ /e " (t '" t)e " iq ' (r " r ' )Tr ^ t ( r '^-^^( r ')^W x ^> +c - c -' 

(153) 

In the case of an experiment in photoluminescence (recombination of photoexcited elec- 
trons and holes) the potential has the form 

V = X)A(r i ,t)-p i , (154) 

3 

and then, after neglecting the photon momentum (dipolar approximation), the luminescence 
spectrum is given by 

P L (v; cu\t)~Yl /k( r > 0/k (r, t)5(Hk 2 /2m x + E G - frw) , (155) 

k 

in arbitrary units, where we have introduced a local approximation and it has been used 
the effective mass approximation for electrons (e) and holes (h), m~ x = m' 1 + m^ 1 is the 
excitonic mass, E G is the energy gap, and f^ h \r,t) are the populations in state k, in 
position r and at time t of electrons (holes) given by 

f: (h \r,t) = [1 + exp{/r(r,t)[ft 2 A; 2 /2m eW - ^(r, t)}' 1 } , (156) 

where /3*(r,t) is the reciprocal of the field of nonequilibrium temperature and fi* e ^(r,t) the 
quasi-chemical potential. 
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VI. ILLUSTRATIVE EXAMPLES 



A. Experiments in ultra-fast-laser spectroscopy. 

Pump-probe experiments in the field of ultrafast laser spectroscopy, devoted to the study 
of the nonequilibrium photoinjected plasma in semiconductors, have been extensively used 
in recent decades, and have been accompanied by a number of theoretical analysis [58-64]. 
In this type of experiments the system is, as a general rule, driven far away from equilibrium 
and consequently, its theoretical description falls into the realm of the thermodynamics of 
irreversible processes in far-from-equilibrium systems, and the accompanying kinetic and sta- 
tistical theories, and a particularly appropriate approach is the nesef described here. The 
theory presented in the previous Section has been used to derive in detail a response func- 
tion theory for the study of ultrafast optical properties in the photoinjected plasma in semi- 
conductors. Particularly, one needs to derive the frequency- and wave number- dependent 
dielectric function in arbitrary nonequilibrium conditions, because it is the quantity which 
contains all the information related to the optical properties of the system (it provides the 
absorption coefficient, the reflectivity coefficient, the Raman scattering cross section, and 
so on). This is described below, and moreover, we describe the application of the results to 
the study of a particular type of experiment, namely the time-resolved reflectivity changes 
in GaAs and other materials [51-54] where signal changes in the reflectivity, AR/R, of the 
order of 1CP 7 are detected, and a distinct oscillation of the signal in real time is observed. 
In Fig. 3 are reproduced time-resolved reflectivity spectra, and in the upper right inset is 
shown the part corresponding to the observed oscillation, as reported by Cho et al. [51] 

Such phenomenon has been attributed to the generation of coherent lattice vibrations, 
and several theoretical approaches have been reported [54-56]. A clear description, on phe- 
nomenological bases, which captures the essential physics of the problem, is reported in Ref. 
[55], and in Ref. [57] is presented an analysis based on nesef, where the different physical 
aspects of the problem are discussed. It is evidenced that the oscillatory effect is provided 
by the displacive excitation of the polar lattice vibrations, arising out of the coupling of 
the carrier-charge density and polar modes, and its decay is mainly governed by the cooling 
down of the carriers. We briefly describe the experiment and all the pertaining NESEF-based 
calculations in continuation. 
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FIG. 3: Reproduction of the time-resolved reflectivity changes in GaAs, as reported by Cho, Kiitt, 
and Kurz in Ref. [51]. 

Let us consider a direct-gap polar semiconductor in a pump-probe experiment. We recall 
that the exciting intense laser pulse produces the so-called highly excited plasma in semi- 
conductors, namely, electron-hole pairs in the metallic side of Mott transition (that is, they 
are itinerant carriers, and we recall that this requires concentrations of these photoinjected 
quasi-particles of order of 10 16 cm~ 3 and up), which compose a two-component Fermi fluid, 
moving in the lattice background. It constitutes a highly nonequilibrated system where the 
photoexcited carriers rapidly redistribute their energy in excess of equilibrium via, mainly, 
the strong long-range Coulomb interaction (pico- to subpico- second scale), followed by the 
transfer of energy to the phonon field (predominantly to the optical phonons, and preferen- 
tially to the LO phonons via Frohlich interaction), and finally via acoustic phonons to the 
external thermal reservoir. Along the process the carrier density diminishes in recombina- 
tion processes (nanosecond time scale) and through ambipolar diffusion out of the active 
volume of the sample (ten-fold picosecond time scale). 

Moreover, a probe interacting weakly with the heps is used to obtain an optical response, 
namely the reflectivity of the incoming laser photons with frequency u and wave vector Q. 
From the theoretical point of view, such measurement is to be analyzed in terms of, as 
already noted, the all important and inevitable use of correlation functions in response 
function theory. The usual application in normal probe experiments performed on a system 
initially in equilibrium had a long history of success, and a practical and elegant treatment 
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is based on the method of the double-time (equilibrium) thermodynamic Green functions 
[45,46]. In the present case of a pump-probe experiment we need to resort to a theory of 
such type but applied to a system whose macroscopic state is in nonequilibrium conditions 
and evolving in time as a result of the dissipative processes that are developing while the 
sample is probed. This is, resorting to the theory of previous section, and we recall that the 
response function theory for nonequilibrium systems needs be coupled to the kinetic theory 
that describes the evolution of the nonequilibrium state of the system. We resort here to 
such theory for the study of the reflectivity experiments of Ref. [51]. 

The time-dependent (because it keeps changing along with the evolution of the macrostate 
of the nonequilibrated system) reflectivity R(ou, Q|t) is related to the index of refraction 
r](cu, Q|t) + ik(uj, Q\t) through the well-known expression 

[,(u,,Q|t) + i] 2 + [«(u,,Qji)] 2 

and the refraction index is related to the time-evolving frequency- and wave vector-dependent 
dielectric function by 

e(u,Q\t) = e'(u,Q\t) + ie"(uj,Q\t) 

= [ri(ou,Q\t) + iK(uj,Q\t)} 2 , (158) 

where rj and e', and k and e", are the real and imaginary parts of the refraction index and 
of the dielectric function respectively. 

We call the attention to the fact that the dielectric function depends on the frequency and 
the wave vector of the radiation involved, and t stands for the time when a measurement is 
performed. Once again we stress that this dependence on time is, of course, the result that 
the macroscopic state of the non-equilibrated plasma is evolving in time as the experiment 
is performed. 

Therefore it is our task to calculate this dielectric function in the nonequilibrium state 
of the system. First, we note that according to Maxwell equations in material media (that 
is, Maxwell equations now averaged over the nonequilibrium statistical ensemble) we have 
that 

£ ->,Q| t )- 1 = I^M (159) 
r(oj, Q) 

where r(ou, Q) is the amplitude of a probe charge density with frequency u and wave vector 
Q, and n(uj, Q\t) is the induced polarization-charge density of carriers and lattice in the 
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media. As shown before the latter can be calculated resorting to the response function theory 
for systems far from equilibrium (the case is quite similar to the calculation of the time- 
resolved Raman scattering cross section [48]), and obtained in terms of the nonequilibrium- 
thermodynamic Green functions, as we proceed to describe. 

Using the formalism we have presented to obtain e(uu, Q\t), [cf. Eqs. (98)] it follows that 

e~\u, Q) - 1 = V(Q) [G cc (ou, Q) + G cl (u, Q) + G ic (uu, Q) + Ga(u>, Q)] , (160) 

giving the reciprocal of the dielectric function in terms of Green functions given by 

(<n c (Q);nt(Q)|a;;f», (161) 

(<n c (Q);nt(Q)|w;t)), (162) 

((n t (Q);nl(Q)\u;t)) , (163) 

«^(Q);nt(Q)|o;;t)>, (164) 

where V(Q) = Aime 2 /VsoQ 2 is the matrix element of the Coulomb potential in plane-wave 
states and n c (Q), and n»(Q), refer to the Q-wave vector Fourier transform of the operators 
for the densities of charge of carriers and the polarization charge of longitudinal optical 
phonons respectively. 

But, the expression we obtain is, as already noticed, depending on the evolving nonequi- 
librium macroscopic state of the system, a fact embedded in the expressions for the time- 
dependent distribution functions of the carrier and phonon states. Therefore, they are to 
be derived within the kinetic theory in nesef, and the first and fundamental step is the 
choice of the set of variables deemed appropriate for the description of the macroscopic 
state of the system. In this case a first set of variables needs be the one composed of the 
carriers' density and energy, and the phonon population functions, together with the set of 
associated Lagrange multipliers that, as we have seen, can be interpreted as a reciprocal 
quasitemperature and quasi-chemical potentials of carriers, and reciprocal quasitempera- 
tures of phonons, one for each mode [42,58,59]. But in the situation we are considering we 
need to add, on the basis of the information provided by the experiment, the amplitudes 
of the LO-lattice vibrations and the carrier charge density; the former because it is clearly 
present in the experimental data (the oscillation in the reflectivity) and the latter because 



G cc (uj, Q) = 

Gdiu;, Q) = 
G ic (u, Q) = 
Ga(u, Q) = 
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of the LO-phonon-plasma coupling clearly present in Raman scattering experiments [60,61]). 
Consequently the chosen basic set of dynamical quantities is 

{H c , N e , N h , n kp , ra£ p , z> q , o q , a q , H B } , (165) 

where 



H, 



C = E[ £ k c kCk + 4^L k ^k] , (166) 

k 

z> q = a q a q , (167) 

Ne = > ^ = X) ^-k ^ k ' ( 16§ ) 

k k 

^k P = 4+pCk , ra£ p = /i_ k _ p /it k , (169) 

with c (c*), ft. (ft*), and a (a*) being as usual annihilation (creation) operators in electron, 
hole, and LO-phonon states respectively (k, p, q run over the Brillouin zone). Moreover, the 
effective mass approximation is used and Coulomb interaction is dealt with in the random 
phase approximation, and then e k = E G + ft 2 |k| 2 /2m e and e£ = ft 2 |k| 2 /2m ft . Finally H B 
is the Hamiltonian of the lattice vibrations different from the LO one. We write for the 
NESEF-nonequilibrium thermodynamic variables associated to the quantities of Eq. (165) 

{&(*), -MWM -0 c (t)iMl(t), F* p (t),F£ p (t), tkuMt), ip^t),^), A,} , (170) 

respectively, where //* and n* h are the quasi-chemical potentials for electrons and for holes; we 
write P c (t) = l//csT*(t) introducing the carriers' quasitemperature T*; /3 q (£) = l/ksT*(t) 
introducing the LO-phonon quasitemperature per mode (o> q is the dispersion relation), 
Po — 1 / IzbTq with T being the temperature of the thermal reservoir. We indicate the cor- 
responding macrovariables, that is, those which define the nonequilibrium thermodynamic 
Gibbs space as 

{E c (t), n{t), n{t), n kp (t), <,(*), i/ q (t), (ajt), (a[\t) = (a q |f)*, E B } , (171) 
which are the statistical average of the quantities of Eq. (165), that is 

E c (t) =Tr{H c g e (t) x g R } , (172) 
n(t) = Tr{N e{h) g e (t) x g R } , (173) 
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and so on, where qr is the stationary statistical distribution of the reservoir and n(t) is the 
carrier density, which is equal for electrons and for holes since they are produced in pairs in 
the intrinsic semiconductor. The volume of the active region of the sample (where the laser 
beam is focused) is taken equal to 1 for simplicity. 

Next step is to derive the equation of evolution for the basic variables that characterize 
the nonequilibrium macroscopic state of the system, and from them the evolution of the 
nonequilibrium thermodynamic variables. This is done according to the generalized ne- 
SEF-based nonlinear quantum transport theory already described, but in the second-order 
approximation in relaxation theory. This is an approximation which retains only two-body 
collisions but with memory being neglected, consisting in the Markovian limit of the theory. 
It is sometimes referred to as the quasi-linear approximation in relaxation theory [43,44], 
a name we avoid because of the misleading word linear which refers to the lowest order in 
dissipation, however the equations are highly nonlinear. 

The NESEF-auxiliary ( "instantaneously frozen" ) statistical operator is in the present case 
given, in terms of the variables of Eq. (165) and the nonequilibrium thermodynamic variables 
of Eq. (171), by 

g(t,0) = exp{-<l>(t)-^(t)[H c -^(t)N e -^(t)N h ]- 

k P 

Y^[^(t)frco^ + ^(t)av + ^(t)a\]-(3 H B } , (174) 
q 

where (pit) ensures the normalization of g(t,0). 

Using such statistical operator the Green functions that define the dielectric function [cf. 
Eq. (160)] can be calculated. This is an arduous task, and in the process it is necessary to 
evaluate the occupation functions 

h(t)=Tr{c k c k g e (t)}, (175) 

which is dependent on the variables of Eq. (171). The (nonequilibrium) carrier quasitem- 
perature T* is obtained, and its evolution in time shown in Fig. 4. 

Finally, in Fig. 5, leaving only as an adjustable parameter the amplitude - which is fixed 
fitting the first maximum -, is shown the calculated modulation effect which is compared 
with the experimental result (we have only placed the positions of maximum and minimum 
amplitude taken from the experimental data, which are indicated by the full squares). 
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FIG. 4: Evolution of the carrier's quasitemperature, calculated in the conditions of the experiment 
in the caption to Fig. 3. 
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FIG. 5: The theoretically evaluated modulation of the time-resolved reflectivity in the conditions 
of Rcf. [51], compared with the experimental data. For simplicity we have drawn only the positions 
of the maxima and minima of the figure in the inset of Fig. 3. 



This demonstrates the reason of the presence of the observed modulating phenomenon 
in the reflectivity spectra, occurring with the frequency of the near zone center LO-phonon 
(more precisely the one of the upper L + hybrid mode [61]) with wave vector Q, the one of 
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the photon in the laser radiation field. The amplitude of the modulation is determined by 
the amplitude of the laser-radiation-driven carrier charge density which is coupled to the 
optical vibration, and then an open parameter in the theory to be fixed by the experimental 
observation. This study has provided, as shown, a good illustration of the full use of the 
formalism of nesef, with an application to a quite interesting experiment and where, we 
recall, the observed signal associated to the modulation is seven orders of magnitude smaller 
than the main signal on which is superimposed. 

B. Charge Transport in Doped Semiconductors 

NESEF is particularly appropriate do describe the transient and steady state of semi- 
conductors in the presence of intermediate to strong electric fields (say tens to hundreds 
of kV/cm), fields which drive the system far- from equilibrium. The question has large 
technological interest because of the presence of such situation in the integrated circuits of 
electronic and optoelectronic devices. 

Let us consider a n-doped direct-gap polar semiconductor, in condition such that the 
extra electrons act as mobile carriers in the conduction band. We use the effective-mass 
approximation, and therefore parabolic band; this implies that in explicit applications it 
needs be controlled the fact that there exists an upper limiting value for the electric field 
strength, such that below this limit we are working with field intensities for which intervalley 
scattering can be neglected. The Hamiltonian of the system is 

H = H + Hi + Han + Hcf + W, (176) 

where 

H = Y.^/^O^W + £ ^( b lK + 1/2) , (177) 

k q,7 

is the Hamiltonian of free electrons and phonons in branches 7 = LO,AC, and 

#i= £ Wq^i+^Ok + M-(q)^4c k+q7 ] , (178) 

k,q,7,o- 

is the interaction Hamiltonian between them. In these equations c(c') and b(b>) are anni- 
hilation (creation) operators in electron states |k >, and of phonons in mode |q 7 > and 
branches 7 =lo,ac (for longitudinal optical and acoustical ones respectively; TO-phonons 
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are ignored once they do not interact with the electrons in the conduction band). Quantity 
M~(q) is the matrix element of the interaction between carriers and 7-type phonons, with 
superscript a indicating the kind of interaction (polar, deformation potential, piezoelectric). 
Moreover, Han stands for the anharmonic interaction in the phonon system, and 

H C F = -Y,eE-Ti, (179) 

i 

is the interaction of the electrons (with charge — e and positions r,) with an electric field F 
of intensity £. The interaction of the system with an external reservoir is taken care of by 
W in Eq. (35); the reservoir is taken as an ideal one - what is satisfactory in most cases - 
and then has its macroscopic (thermodynamic) state characterized by a canonical statistical 
distribution with temperature T . 

Consider now the nonequilibrium thermodynamic state of the system: the presence of 
the electric field changes the energy of the electrons (they acquire energy in excess of equi- 
librium), and these carriers keep transferring this excess to the lattice and from the lattice 
to the thermal reservoir, and an electrical current (flux of electrons) follows. Thus, we need 
to choose as basic variables 

{E e (t), N e (t), P e (t), E LO (t), E AC (t),E R }, (180) 

that is, respectively, the energy, number, and linear momentum of the carriers, the energies 
of the LO and AC phonons, and the energy of the reservoir; the latter is constant in time for 
being considered as an ideal one. The corresponding dynamical quantities are 

{H e ,N e ,P e ,H LO ,H AC ,H R }, (181) 

i.e. the Hermitian operators for the partial Hamiltonians, the electron number and the 
linear momentum. We noticed that the above choice implies in disregarding electro-thermal 
effects, whose inclusion would require to introduce the flux of energy (heat current) of the 
carriers; it has a minor influence on the results to be reported. 

According to the nonequilibrium statistical ensemble formalism described in section 2 of 
the preceding article, the nonequilibrium thermodynamic state of the system, in an alter- 
native description to the one provided by the variables of Eq. (181), can be completely 
characterized by a set of intensive nonequilibrium thermodynamic variables (Lagrange mul- 
tipliers that the variational construction of the formalism provides), namely 

{F e (t), F ne (t), F e (t), F LO (t), F AC (t), f3 }}, (182) 
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The variables in this Eq. (182) are present in the auxiliary statistical operator which the 
formalism introduces, in this case given by 

g(t, 0) = exp{-0(t) - F ne (t)N e - F e {t)H e - 

-F e (i) • P e - F LO (t)H LO - F AC (t)H AC } x g R }, (183) 

where g B is the canonical distribution of the reservoir at temperature T . We recall that the 
operator of Eq. (183) is not the statistical operator describing the macroscopic state of the 
system, which is a superoperator of this one, and <fi(t) (playing the role of a logarithm of a 
nonequilibrium partition function) ensures the normalization of g(t,0) (see Appendix A). 
The intensive nonequilibrium thermodynamic variables of Eq. (182) are usually redefined 

as 

F e (t) = P* e (t) = [k B T:(t)]-\ (184) 
F ne (t) = -P* e (t)ti* e (t), (185) 

F e (t) = -&(t)v e (t), (186) 

F LO (t) = P* LO (t) = [k B T* LO (t)]~\ (187) 

F AC (t) = /3* AC (t) = [k B T* AC (t)}-\ (188) 

and we recall that /3 [in Eq. (182)] is [faTo] -1 . These Eqs. (184) to (188) introduce 
the so-called quasitemperatures, T*(t), T£ (t), T AC (t), of electrons and phonons, and the 
quasi-chemical potential fJ>l(t) and the drift velocity v e (i) of the electrons; k B is as usual 
Boltzmann constant. 

Proceeding with the calculations of the equations of evolution of the basic variables, the 
corresponding set of equations of evolution are obtained, which have expressions of the form 
[42,59,62] 

j t E e (t) = -^E.P e (t)-J^(t), (189) 

jP e (t) = -nVeE + jg(t) + J { Zm P (t) , (190) 
f t E LO (t) = 4l(t) - J^ AN (t) , (191) 

jE AC (t) = jf AC (t) + J^ AN (t) - J A \ dif {t) , (192) 

where, we recall, E e (t) is the carriers' energy and P e (t) the linear momentum; ELoif) the 
energy of the LO phonons which strongly interact with the carriers via Frohlich potential 
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in these strong-polar semiconductors (hence it predominates over the nonpolar-deformation 
potential interaction and then the latter is disregarded); £Uc(0 is the energy of the acoustic 
phonons playing the role of a thermal bath; and E stands for the constant electric field. 

Let us analyze these equations term by term. In Eq. (189) the first term on the right 
accounts for the rate of energy transferred from the electric field to the carriers, and the 
second term accounts for the transfer of the resulting excess energy of the carriers to the 
phonons. In Eq. (190) the first term on the right is the driving force generated by the 
presence of the electric field. The second term is the rate of momentum transfer due to 
the interaction with the phonons, and the last one is a result of scattering by impurities 
(these two terms are then momentum relaxation contributions). In Eq. (191) and Eq. (192) 
the first term on the right describes the rate of change of the energy of the phonons due to 
interaction with electrons. More precisely they account for the gain of the energy transferred 
to them from the hot carriers and then the sum of contributions Je\ q (t) an d Je\ c (*) * s ec L ua l 
to the last in Eq. (189), but accompanied with a change of sign. The second term in Eq. 
(191) accounts for the rate of transfer of energy from the optical phonons to the acoustic 

(2) 

ones, via anharmonic interaction. The contribution J L q AN (t) is the same but with different 
sign in Eq. (191) and Eq. (192). Finally, the diffusion of heat from the AC phonons to the 
reservoir is accounted for in the last term in Eq. (192). The detailed expression for the 
collision operators are given in Ref. [42] (quantities J are positive). 

The solution of these equations allows for a detailed analysis of the nonequilibrium ther- 
modynamic state and transport properties of these materials. Let us consider first the 
Ill-Nitride compounds, which nowadays present a particular interest as a result of their po- 
tential use in lasers and diodes emitting in the blue and ultraviolet region (see for example 
[63]). 

Let us consider the steady state which follows very rapidly (in a hundred-fold femtosec- 
ond time scale), what can be understood on the basis of the action of the intense Frohlich 
interaction in these strong polar semiconductors, with the rate of transfer of energy from 
carriers to LO phonons rapidly equalizing the rate of energy pumping from the external field 
of intensity £, even at high fields. It is then characterized by the constant-in-time variables 
quasitemperature, T*, drift velocity, v e , and quasi-chemical potential, //*, all referring to 
the electron system, and T£ G , the quasitemperature of the LO phonons, and ta, the qua- 
sitemperature of the acoustic phonons. In Fig. 6 it is shown the dependence of T* with the 
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electric field, while in Fig. 7 is presented such dependence for the drift velocity for n-doped 
GaN, with n = 10 17 cm~ 3 . The quasi-chemical potential is determined by the values of the 
concentration and the electron quasitemperature, and the deviation of T£ and T4 from the 
value in equilibrium is small and can be neglected. 
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FIG. 6: Dependence of the quasitemperature T* with the electric field in the steady state of n- 
doped GaN. We have used m* = 0.19mo, £0 = 9.5, a carrier concentration n = 10 17 cm~ 3 , and the 
thermal bath temperature is To = 300 K. After Rcf. [64] 

We can now proceed to calculate the cross section for scattering of light by electrons in 
the electric field S. Taking into account that the interaction of electrons and radiation is 
given by 

VerQ, = G(Q)n Q + H.c. , (193) 

where HQ is the momentum of the photon, G(Q) is the matrix element of the interaction, 
and 

™Q = Yl ™ k Q = Yl C ic+Q c k , (194) 

k k 

the cross section of Eq. (143), in first order in the scattering operator, is given by 



d 2 a(Q,u) V 2 fko 2 



dudtt (2ttH) 3 c 4 



\G(Q)\ 2 / dre-^Tr{n' Q n Q n £ } , 



(195) 



and we recall that the system is in a steady state, and, moreover, both integrals in Eq. (143) 
can be combined in the given form above, and we have used that g(cu) = h 3 u 2 /c 3 and vq = c. 
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ELECTRIC FIELD (kV/cm) 

FIG. 7: Dependence of the drift velocity of the electrons with the electric field in the steady state 
of n-doped GaN, in the same experimental conditions indicated in the caption to Fig. 6. After 
Ref. [64]. 

We can now use the generalized fluctuation-dissipation theorem of the preceding section. 
Using the approximation of keeping only what we have termed as the relevant part, i.e., 
taking TZ of Eq. (183) instead of 1Z S) in the steady state, we find that 

f^S^) e -0(**-Q-Ve)]-l /me -l(Q >w |£) ) (196) 

dudil 

In this Eq. (196), e(Q, u)\£) is the dielectric function at wave vector Q and frequency u (the 
momentum and energy transfer in the scattering event as we have seen), once we use that 

where in the effective mass approximation £(k) = h 2 k 2 /2m*, and Eq. (197) is of the form 
of Lindhardt (RPA) dielectric function, but in terms of the nonequilibrium distribution 
functions f(k\£). They have an expression of a drifted Fermi-Dirac-like distribution (with 
the presence of the electric field-dependent quasitemperature and quasi-chemical potential), 
which in the usual experimental condition can be appropriately approximated by a drifted 
Maxwell-Boltzmann-like distribution, namely 

f(k\£) = AeM-mim) — E ■ Kk] } , (198) 

with 

«< f) =wW (199) 
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A{£) 



87i 3 nh 3 [f3:{£)]V 2 



(200) 



(27rm*) 3 / 2 

where it has been used that V(Q) = Aire 2 /(Ve Q 2 ) is the Fourier transform of the Coulomb 
potential with e being the static dielectric constant and V the volume of the system. More- 
over, s is a positive infinitesimal which is taken in the limit of going to +0 to produce the 
real and imaginary parts of e(Q,o>) = ei(Q, u>) + ie 2 (Q,uj). Going over the continuum, i.e., 
transforming the summation in Eq. (197) in an integral and using spherical coordinates k, 
9, ip we find for the real part of the dielectric function [65] 

, V(Q) 



ei(Q,w) 



where 



with 



F(Q,u>) 



47T 3 

nn 3 y/2 5 m* e l3* 

hQ 



[D( Vl ) + D(y 2 )} , 



D(y) 

which is Dawson's integral, and 



y 

exp(— y 2 ) J exp(x 2 )dx . 
o 



(201) 



(202) 



(203) 
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2 + ng ( "- Q ' Ve) 



Q m* e 
~2~frQ 



[u-Q- v e ) 



(204) 



(205) 



Substituting Eq. (202) in Eq. (201), we obtain that 



where 



^i(Q^) = i + i/^^[^(yi) + ^)], 

2 47re 2 n 



(206) 



(207) 



is the Debye-Huckel screening factor. On the other hand, the imaginary part is given by 

(208) 



e 2 (Q^) = n f^k H \zM-yl) - exp(-^ 2 )] . 



2f3* e hQ 3 

Let us consider the case of GaN in the steady-state thermodynamic conditions indicated 
in Figs. 6 and 7, in the case: £ = 100 kV/cm, and three experimental geometries, i.e. three 
values of the scattering angle 9. The scattering spectrum is composed of two contributions, 
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FIG. 8: The Raman band due to scattering by the single-particle elementary excitations: £ = 100 
kV/cm, Q = 1.8 x 10 5 cm _1 , the thermodynamic state as characterized by Figs. 6 and 7 for the 
given £ , and the three scattering angles, 6 = (full line), 6 = ir/4 (dashed line), and 6 = it/ 2 
(dotted line). After Ref. [65]. 

namely, scattering, in Fig. 8 with Q — 1.8 x 10 5 cm -1 , by individual electrons at low ou, and 
in Fig. 9 scattering, with Q = 5 x 10 4 cm -1 , by collective excitations (plasmons) around 
the plasma frequency u? v = Airne 2 /e m* e . In these Figs. 8 and 9 it can be noticed a shift 
in the scattering bands, depending on the scattering angle. This is a result of the presence 
of the term Q • v e in Eq. (196), which is null for the experimental geometry in which the 
momentum transfer Q is perpendicular to the drift velocity, which is in the direction of 
the electric field (v e = .M e E, where Ai e is the electron mobility), and maximum when Q 
and v e are parallel. This has a quite interesting consequence, consisting in that the shift 
in frequency permits a measurement of the drift velocity, and then of the mobility in the 
conditions of the experiment. If we call Acj pea k the difference of frequencies at the peak 
positions of the bands for scattering by plasmons, for 6 = o and 9 = n/2 (Fig. 9), then 
v e = Aup ea k/Q and M e = v e /S = Au peak /QS. 
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FIG. 9: The Raman band due to scattering by the plasma excitation: same specifications as in the 
caption to Fig. 8, but for Q = 5 x 10 4 cm _1 . After Ref. [65]. 

VII. FINAL REMARKS 

We have presented a Response Function Theory, accompanied of a Fluctuation Dissipa- 
tion Theorem, and a Theory of Scattering adapted to deal with systems arbitrarily away 
from equilibrium, including situations of time and space experimental resolution. 

Such theory was built within the framework of a Gibbs-style Non-Equilibrium Statistical 
Formalism. The general form of the generalized susceptibility (space and time dependent) 
is obtained in the form of space and time dependent correlation functions defined over the 
nonequilibrum ensemble. Moreover it is dependent on the variables that characterized the 
nonequilibrium thermodynamic state of the system. Therefore, the generalized suscepti- 
bility of the corresponding experimental situation is coupled to the equations of evolution 
of the nonequilibrium thermodynamic variables. It is also presented the method eventu- 
ally useful for calculations, of nonequilibrium thermodynamic Green functions, that is, the 
extension to arbitrary nonequilibrium conditions of Bogoliubov-Tyblikov thermodynamic 
Green functions. 

A Fluctuation-Dissipation theorem in the context of the Nonequilibrium Ensemble For- 
malism has been presented in Section III (see also subsection IV. A), which is an extension 
to arbitrary nonequilibrium conditions of Kubo's one. 
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In Section V we have presented a Theory of Scattering appropriate for scattering ex- 
periments done on systems which are arbitrarily away from equilibrium. A space and time 
dependent scattering cross section is obtained, which, as in the case of the generalized sus- 
ceptibility of Response Function, is dependent on the nonequilibrium thermodynamic state 
of the system, and then the scattering cross section is coupled to the equations of evolution 
of the nonequilibrium variables. 

Finally, in section VI we have presented some illustrative examples of the working of the 
theory in the analysis of several experimental situations, namely, in ultrafast laser spec- 
troscopy and charge transport in doped polar semiconductors. 

Acknowledgments: The authors would like to acknowledge partial financial support 
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Appendix A: The Nonequilibrium Statistical Operator 

Construction of nonequilibrium statistical ensembles, that is, a Nonequilibrium Statistical 
Ensemble Formalism, NESEF for short [8-11], consisting in, basically, the derivation of a 
nonequilibrium statistical operator (probability distribution in the classical case) has been 
attempted along several lines. In a brief summarized way we descrive the construction of 
NESEF within a heuristic approach, and, first, it needs to be noticed that for systems away 
form equilibrium, several important points need to be carefully taken into account in each 
case under consideration: 

1. The choice of the basic variables (a wholly different choice than in equilibrium when 
it suffices to take a set of those which are constants of motion), which is to be based 
on an analysis of what sort of macroscopic measurements and processes are actually 
possible, and moreover, one is to focus attention not only on what can be observed 
but also on the character and expectative concerning the equations of evolution for 
these variables [11,66]. We also notice that even though at the initial stage we would 
need to introduce all the observables of the system, an eventually variances, as time 
elapses more and more contracted descriptions can be used when it enters into play 
Bogoliubov's principle of correlation weakening and the accompanying hierarchy of 
relaxation times [67]. 

2. The question of irreversibility (or Eddington's arrow of time) on what Rudolf Peierls 
stated that: "In any theoretical treatment of transport problems, it is important 
to realize at what point the irreversibility has been incorporated. If it has not been 
incorporated, the treatment is wrong. A description of the situation that preserves the 
reversibility in time is bound to give the answer zero or infinity for any conductivity. 
If we do not see clearly where the irreversibility is introduced, we do not clearly 
understand what we are doing" [68]. 

3. Historicity needs be introduced, that is, the idea that it must be incorporated all the 
past dynamics of the system (or historicity effects), all along the time interval going 
from a starting description of the macro-state of the sample in the given experiment, 
say at to, up to the time t when a measurement is performed. This is a quite important 
point in the case of dissipative systems as emphasized among others by John Kirkwood, 
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Green, Robert Zwanzig and Hazime More [15-19]. It implies in that the history of the 
system is not merely the series of events in which the system has been involved, but it is 
the series of transformations along time by which the system progressively comes into 
being at time t (when a measurement is performed), through the evolution governed 
by laws of mechanics. [69] 

Concerning the question of the choice of the basic variables, differently to the case in equi- 
librium, immediately after the open system of N particles, in contact with external sources 
and reservoirs, has been driven out of equilibrium, it would be necessary to describe its state 
in terms of all its observables and, eventually, introducing direct and cross-correlation. But, 
as time elapses Bogoliubov's principle of correlation weakening allow us to introduced in- 
creasing contractions of descriptions. Let us say that we can introduce a description based on 
the observables {Pj}, j = 1, 2, n, on which depends the noneuilibrium statistical operator. 

On the question of irreversibility Nicolai S. Krylov [21] considered that there always exists 
a physical interaction between the measured system and the external world that is constantly 
"jolting" the system out of its exact microstate. Thus, the instability of trajectories and the 
unavoidable finite interaction with the outside would guarantee the working of a "crudely 
prepared" macroscopic description. In the absence of a proper way to introduce such effect, 
one needs to resort to the interventionist's approach, which is grounded on the basis of such 
ineluctable process of randomization leading to the asymmetric evolution of the macro-state. 

The "intervention" consists into introducing in the Liouville equation of the statistical op- 
erator, of the otherwise isolated system, a particular source accounting for Krylov's "jolting" 
effect, in the form (written for the logarithm of the statistical operator) 

^ In M £ (t) + 1 [In » e (i), H] = -epn tt e (i) - In 0)], (Al) 

where e (kind of reciprocal of a relaxation time) is taken to go to +0 after the calcula- 
tions of average values has been performed. Such mathematically inhomogeneous term, in 
the otherwise normal Liouville equation, implies in a continuous tendency of relaxation of 
the statistical operators towards a referential distribution, $t, which, as discussed below, 
represents an instantaneous quasi-equilibrium condition. 

We can see that Eq. (A.l) consists of a regular Liouville equation but with an infinitesimal 
source, which provides Bogoliubov's symmetry breaking of time reversal and is responsible 
for disregarding the advanced solutions [8,11,70]. This is described by a Poisson distribution 
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and the result at time t is obtained by averaging over all t' in the interval (t ,t), once the 
solution of Eq. (A.l) is 



tt e (f) = exp | -S(t, 0) + J dt'e'V-O-^Stf, t' - t)\ , (A2) 



to 



where 



S(t,0) = -\n^(t,0), (A3) 

S(t', t'-t) = exp | - ^r(*' - t)#i S(t', 0) exp (^(f' - t)H[\ , (A4) 
( w J J 

and the initial-time condition at time to, when the formalism begins to be applied, is 

n e (t ) = 0). (A5) 

In and 5, the first time variable in the argument refers to the evolution of the nonequi- 
librium thermodynamic variables and the second to the time evolution of the dynamical 
variables, both of which have an effect on the operator. 

This time t , of initiation of the statistical description, is usually taken in the remote 
past (t — > — oo) introducing an adiabatic switching-on of the relaxation process, and in Eq. 
(A. 2) the integration in time in the interval (to,t) is weighted by the kernel exp{e(t f — t)}. 
The presence of this kernel introduces a kind of evanescent history as the system macro-state 
evolves toward the future from the boundary condition of Eq. (A. 5) at time (£ — > — oo) 
a fact evidenced in the resulting kinetic theory [8-11,14,17] which clearly indicates that it 
has been introduced a fading memory of the dynamical process. It can be noticed that the 
statistical operator can be write in the form 

ft e (f) = &(i,0) ■ ( A6 ) 

involving the auxiliary probability distribution §t(t, 0), plus $H e (t) which contains the his- 
toricity and irreversibility effects. Moreover, in most cases we can consider the system as 
composed of the system of interest (on which we are performing an experiment) in contact 
with ideal reservoirs. Thus, we can write 

&(t, 0) = p(t, 0)x PR . (A7) 
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and 

9ft e (t) = p £ {t) x p R , (A8) 

where p e (t) is the statistical operator of the nonequilibrium system, p the auxiliary one, and 
Pn the stationary one of the ideal reservoirs, with p £ (t) given then by 

Pe(t) = exp j-S(f,0) + J dt'e^'-V^Sit'J j , (A9) 

having the initial value p(t ,0) (to — >■ _ oo), and where 

S{t,0) = -mp(i,0) , (A10) 

Finally, it needs be provided the auxiliary statistical operator p(t,0). It defines an in- 
stantaneous distribution at time t, which describes a "frozen" equilibrium defining at such 
given time the macroscopic state of the system, and for that reason is sometimes dubbed 
as the quasi- equilibrium statistical operator. On the basis of this (or, alternatively, via 
the extremum principle procedure [11,71-76], and considering the description of the non- 
equilibrium state of the system in terms of the basic set of dynamical variables Pj, the 
reference or instantaneous quasi-equilibrium statistical operator is taken as a canonical-like 
one given by 

n 

p(t, 0) = exp{-0(t) - F j(t) p j} , (AH) 

j 

in the classical case, with <j)(t) ensuring the normalization of p and playing the role of a kind 
of a logarithm of a partition function, say, <p(t) = \nZ(t). Moreover, in this Eq. (A. 11), Fj 
are the nonequilibrium thermodynamic variables associated to each kind of basic dynamical 
variables Pj. The nonequilibrium thermodynamic space of states is composed by the basic 
variables {Qj(t)} consisting of the averages of the {Pj} over the nonequilibrium ensemble, 
namely 

Q j (t)=Tr{P jPe (t)}, (A12) 
which are then functional of the {Fj(t)} and there follow the equations of state 

where 5 stands for functional derivative. 
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Moreover 

S(t) = Tr{S(t, 0)p(t, 0)} = -Tr{p(t, 0) lnp(t, 0)} , (A14) 

is the so-called informational entropy characteristic of the distribution p, a functional of the 
basic variables {Qj(t)}, and it is verified the alternative form of the equations of state given 

- SI - w ■ < ais > 
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